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We show that the classical equations of gravity follow from a thermodynamic rela- 
^ I tion, 6Q = T6S, where S is taken to be the Wald entropy, applied to a local Rindler 

■ horizon at any point in spacetime. Our approach works for all diffeomorphism- 

. invariant theories of gravity. This suggests that classical gravity may be thermody- 



namic m origm. 



INTRODUCTION 



r-| I Black holes have often provided insights into the nature of quantum gravity and the 

structure of spacetime. For example, the holographic principle emerged through Gedanken- 
experiments that took a property of black holes — the sub-extensive scaling of Bekenstein- 
^ ■ Hawking entropy — and applied it to arbitrary gravitational systems. 

Another potentially profound insight comes from a paper by Jacobson Jacobson 

m 

O ' considered the puzzling fact that the laws of black hole mechanics, derived in classical general 
O ' relativity, seem mysteriously to anticipate the laws of black hole thermodynamics, derived in 
- ^ , semi-classical gravity. Rather than trying to explain how classical laws could "know about" 
^ ■ quantum-mechanical ones, Jacobson reversed the logic, regarding the thermodynamics to be 
a premise rather than a consequence. Quite remarkably, by assigning the thermodynamic 
properties of black hole horizons to local light-cones in spacetime (not necessarily near a 
black hole), the Einstein equation re-appears as an equation of state. This seems to suggest, 
as indeed many theorists believe, that gravity is not a fundamental theory but originates in 
some kind of thermodynamic approximation. 

The question arises whether this alluring result is somehow an artifact of Einstein grav- 
ity, or whether the connection between thermodynamics and gravity goes deeper, persisting 
also in general, higher-curvature theories of gravity. But extending the original derivation 
to higher-curvature theories is nontrivial, in part because that derivation makes use of the 
Raychaudhuri equation, whose usefulness is obscured in higher-curvature theories: the Ray- 
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chaudhuri equation relates the derivative of the expansion of the horizon to the Ricci tensor, 
but a simple relation between the Ricci tensor and the stress tensor holds only for Einstein 
gravity. Moreover, in generic theories of gravity, the entropy is not simply proportional to 
the area. 

In this paper, we obtain the classical gravitational equations from thermodynamics with- 
out making use of the Raychaudhuri equation. Specifically, we show that the classical 
equations of gravity follow directly from the Clausius relation, 6S = 6Q/T. Here for S we 
use Wald's definition of entropy, which is the entropy (in place of A/ 4) that satisfies the 
first law of thermodynamics in higher-curvature theories. Our result suggests that classical 
gravity might have a quite intriguing thermodynamic origin. 



GENERAL THEORIES OF GRAVITY 



Consider a general diffeomorphism-invariant theory of gravity in any number of dimen- 
sions. For convenience, we will assume that the Lagrangian is a polynomial in the Riemann 
tensor but does not involve its derivatives. One may regard the Lagrangian formally as 
dependent on both the metric and the Riemann tensor even though of course the Riemann 
tensor, depends on the metric . Specifically, let the action be 

1 f 



d X^/^L{gab: Rabcd) ■ (l) 
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We have set Newton's constant to unity. Define 

pabcd _ '-'^ 

dRabcd 

pabcd j^g^g ^YiQ same algebraic symmetries as the Riemann tensor, including cyclicity. One 
then finds that the equation of motion that follows from ([T]) (supplemented by appropriate 
generalizations of Gibbons-Hawking-like boundary terms) is 

Pa'^'Rbcde - ^V'VPacdb ' ^Lgab = MT^b ■ (3) 

For example, when the Lagrangian is L = /(-R), we find P'^^'^ = ^f'{R) (^g^^g^<^ — g"-^g^^Y 
Thus, the equation of motion is 

f'{R)Rab - VaV J'(i?) + (af'{R) - IfiR)] gab = SnT^b . (4) 
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This reduces to Einstein's equation when fiR) = R- 

Another example is Lovelock gravity 4, |5|, the most general extension of Einstein gravity 
for which the equations of motion do not contain derivatives of the Riemann tensor. The 
Lagrangian is L = J2m=Gi ^mLm^ where Cm are constants of dimension (length) which 
are arbitrary as far as gravity is concerned, and mmax = (-D — 2)/2 for even D dimensions 
and m^ax = {D — l)/2 for odd D. Each term is made up of contractions of products of 
the Riemann tensor: 

r _ ^ S:il---i2m pili2 ■phm-\32m f t\ 

2m h-j2m^ ni2---^ i2m-li2m " > 

Here the 5 symbol is the generalized Kronecker delta, defined as the sum over signed per- 
mutations of products of ordinary Kronecker deltas. The Einstein-Hilbert action with a 
cosmological constant is just a special case of the Lovelock action with ci = 1 and cq = — 2A. 
When D < 4, there are no other possible terms; the next term appears for D > 5. It is 
L2 = R^ — 'iR°'^Rab + R°'^"^Rabcd^ known as Gauss-Bonnet gravity, which appears in the 
low-energy effective action of certain string theories [g, 7]; its coefficient in ten-dimensional 
heterotic string theory is C2 = +a' /A. The Gauss-Bonnet action is a topological invari- 
ant in four dimensions, just as the Einstein-Hilbert action is a topological invariant in two 



dimensions. It is convenient to write ([5]) in the form js] 



Lm — Qlm) Rabcd ■ (6) 



Then P'^*'^'^ = mQ"^p which has the nice property that VaP"'^^'^ = 0. The equation of 
motion for Lovelock theory is 

^ c„ ( m Q^^^i?",,, - -L^g^' ] = SttT'^" , (7) 

m=0 ^ ^ 

which follows easily from ([3]). 

In each of these theories, one can associate an entropy with Killing or black hole horizons. 
For example, in place of A/ 4, the entropy in f{R) gravity is 

Sf = fiR) J , (8) 
while for Gauss-Bonnet gravity, black holes have an entropy of 



i J rf^-2a;v^(l + 2c2(^-^)i?) , (9) 
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where '^^R is the scalar curvature of (the cross-section of) the horizon. We will show below 

n 

that, as in Jacobson's derivation of Einstein's equation from 5* = A/ 4 [1], varying these 
entropies and imposing the Clausius relation, 6Q = T6S, leads directly to the equations of 
classical gravity. 



WALD ENTROPY 



Wald 0] and other authors [ol, Q have developed a powerful and elegant Lagrangian- 
based method for determining the entropy of a black hole with a Killing horizon. Wald's 
method works for any diffeomorphism-invariant theory in any number of dimensions and 



does not require Euclideanization. Here we adopt a simplified version of the formalism [11 1. 
Consider a generally covariant Lagrangian, L, that depends on the Riemann tensor but does 
not contain derivatives of the Riemann tensor. Under the diffeomorphism x"" ^ x"" + the 
metric changes via 6gab = —^a^b — ^bC,a- By diffeomorphism-invariance, the change in the 
action, when evaluated on-shell, is given only by a surface term. This leads to a conservation 
law, VqJ" = 0, for which we can write J"^ = VbJ""^, where J°'^ defines (not uniquely) the 
antisymmetric Noether potential associated with the diffeomorphism 2|]. 

For a Lagrangian of the type L = L{gab, Rabcd) direct computation shows that J"-^ is given 
by (see |ll[ |) 



jab ^ _2pabcdy^^^ ^ (V.P"''"'^) , (10) 
with P"-^'^''- = dL/dRabcd- The Noether charge associated with a rigid diffeomorphism is 



defined by integrating the Noether potential over a closed spacelike surface S: 



Q 



J"'^dSab 



(11) 



When is a timelike Killing vector (the one whose norm vanishes at the Killing horizon) , it 
turns out [2|, yj that the corresponding Noether charge is precisely the entropy, S, associated 



with the horizon, apart from a few factors: 



S 



dSahJ"'^ 



(12) 



Here k is the surface gravity of the black hole horizon. The integral for this "Wald entropy" 
can be evaluated over any spacelike cross-section of the Killing horizon 9|]. In fact we can 
formally define the quantity S on any closed spacelike surface, S, of codimension two (such 
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as a section of a stretched horizon), and only at the end take the limit in which that S 
approaches a section of the Killing horizon. It can be shown, for example, that both ([8]) and 
([9]) are just special cases of Wald entropy. 



GRAVITATION FROM THERMODYNAMICS 



Now let us show how the classical equations of gravity, ([3]), arise thermodynamically. 
(That the equations look thermodynamical has been shown for spherically-symmetric Love- 
lock gravity [l2^.) The set-up is as follows Take any spacetime point p and pick any 
future-directed null vector fc" emanating from p. In the vicinity of p, the plane orthogonal to 
k"" defines a local acceleration, or Rindler, horizon, H. Let Bi be any spacelike neighborhood 
of p of codimension two that locally lives on the Rindler plane, and let B2 be some further 
section of the Rindler plane along k"-. Next, let be a future-directed approximate timelike 
Killing vector that generates boosts and asymptotically approaches k"". The orbits of and 
the plane orthogonal to the acceleration vector of C,"' define a stretched horizon, S. As in 
the membrane paradigm 1^, [l^, points on H and points on S can be put in one-to-one 
correspondence by, say, ingoing null rays that pierce both surfaces. Let Si be the images of 
Bi on S via this correspondence. See Fig. 1. 

Let ^"^a = — a^, where the norm a (which turns out to be a lapse) is taken to be constant 
over S. This norm vanishes at H, a Killing horizon. Let u"' be the proper velocity of a fiducial 
observer moving along the orbit of i.e. u'^ = (^)° = where r is the proper time. Let 
n°' be the spacelike unit normal to S, pointing in the direction of increasing a. Both and 

map to k"- in the limit that a ^ 0, for which H. 

After these preliminaries, we are ready to deduce the classical equations of gravity from 
thermodynamics. The key idea {l| is to assign black hole thermodynamic properties to local 



Rindler horizons 



151]. The stretched horizon is assigned a local temperature, Tioc = K/27ra, 



as well as the Wald entropy appropriate to the given theory of gravity. 

By ( [T0|) and (fT2|) . the Wald entropy associated with a stretched horizon at time r is 



S = -^ [ dSat {P^'^'V,^, - 2^,V,P^'''') . (13) 
Js{t) 
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FIG. 1: The local Rindler horizon, i7, of an arbitrary spacetime point p is defined by a null vector, 
A;". A stretched horizon, S, is defined by a timelike approximate Killing vector ^" and has a normal 

vector field n". Bi and Si are spacclikc patches of codimcnsion two that inhabit the planes of the 
Rindler and stretched horizons in the directions orthogonal to the figure, with p contained in B\. 

Next, we vary the entropy along the timelike congruence. The entropy change is 



5S = 5(t2)-5(ti) 

1 

4k 



JS{t2) JS{ti) 

= d-LaVb {P^'^'VcU - 2CdVeP"'^') (14) 

In the last step, we have used Stokes' theorem for an antisymmetric tensor field A^^: 



dEaVbA''" = - dSabA''" , (15) 
E Js 

where our E has the boundary S = S{ti) U S{t2), and the minus sign comes about because 
E is timelike. (To be explicit, our conventions here are dT^a — dAdr and dSab — ^{naUb — 
Uanb)dA, where the normal n"' to the stretched horizon points outwards, away from the 
true horizon.) Recall that P"'^'^ has the same algebraic symmetries as the Riemann tensor, 
including cyclicity. Using those symmetries, we find that 

So far we have not used any properties of Now we will assume that ^" is an approximate 
timelike KiUing vector. An exact KiUing vector satisfies KiUing's equation V^Cc + VcCb = 
from which it follows that VaV^^c = R'^abc^d- An approximate Killing vector indeed satisfies 
Killing's equation locally. We will also assume the apphcability of the second equation 
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within our local Rindler patch, a point we will discuss in the next section. The terms in 
parentheses drop out by symmetry. Using our assumption, we find 

TlocSS =^ [ {P'^'^^RdcbeC- 2UVbVcP'''"') UadrdA (17) 
Sira is 

On the other hand, the locally-measured energy or heat fiux into the stretched horizon is 

6Q = + [ (HlaT^u^ = - [ dAdrUaT^C ■ (18) 

Now S is not yet the entropy of the true horizon, H, since we still have to take the limit 
in which the stretched horizon becomes null. Then both n'^ and become proportional to 
the null vector k"- (with the same proportionality constant). Equating 6Q and Ti^cSS and 
taking the null limit, we obtain 

Since this holds for all null vectors fc" at p, we infer that 

Pa'^'Rbcde - 2W^Pacdb + ^ 9ab = STcTab , (20) 

for some scalar function (f. By demanding conservation of the stress tensor and using the 
Bianchi identities, we find that ip = — |L + A, where A is an integration constant. Thus we 
see that imposing TiocSS = 6Q at any point in spacetime necessarily implies that 

P -fi., - 2V'V^P^. - + Ag^ = S.T^ . (21) 

With the cosmo logical constant appearing as an integration constant, this is precisely the 
classical equation of motion, ([3]), for our theory of gravity. 

DISCUSSION 

We have shown that the equations of classical gravity follow from thermodynamics. Our 
derivation did not require the Raychaudhuri equation. Moreover, since we started with 
the Wald entropy, we could go beyond the Einstein equation to the equations of motion of 
general theories of gravity. Since these were obtained from the Clausius relation, they can 
be regarded as equations of state — relations between thermodynamic state variables. 
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Satisfying as this is, there remain some loose ends. One observation [IGy is that the 
derivation of the Wald entropy itself relies on the equations of motion being obeyed. Al- 
though our approach never explicitly invokes the equations of motion, it is still unclear 
whether any derivation, including Jacobson's original calculation, that begins with an on- 
shell expression which agrees with Wald entropy (such as A/ 4) is implicitly assuming the 
answer, or whether that is simply self-consistency. Another technical concern is our use of 
the equation VaV^^c = Rabcd^'^- This equation is obviously true when C,"' is an exact Killing 
vector but to what extent can one trust it when is an approximate Killing vector? For a 
general spacetime, Riemann normal coordinates can be applied to any local patch. In such 
coordinates an approximate boost Killing vector looks like xdt + tdx- However, such a vector 
does not obey VaV^^c = Rabcd^'^ in general, so our assumption cannot be satisfied through 
coordinate choices alone. In this light, it is interesting that for f{R) theories, previous work 
has found that the Clausius relation does not give the equations of motion but also has 



additional terms 17|]; these have been interpreted as non-equilibrium effects. Perhaps the 
failure of VaVb^c = RabcdC,'^ to hold may be traced to such effects. In that case, our deriva- 
tion may indicate how to determine potential non-equilibrium terms for a general theory of 
gravity. (On the other hand, we also make use of this equation for Einstein gravity where 
there are no such terms, so perhaps it is an innocuous assumption.) It would be interesting 
to understand this better, as well as to connect our method to previous approaches [itI . 



While this paper was being prepared, the preprint 181] appeared, claiming similar con- 
clusions; unfortunately, among other things, their starting formula for entropy (equation 9 
in [3]) is manifestly incorrect, leaving the result in doubt. 
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